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I I Abstract: Deep-inelastic production of heavy quarks at HERA, especially charm, 

is an excellent signal to measure the gluon distribution in the proton at small x 
' values. By measuring various differential distributions of the heavy quarks this 

>• ■ reaction permits additional more incisive QCD analyses due to the many scales 

^ , present. Furthermore, the relatively small mass of the charm quark, compared to 

^ I the typical momentum transfer Q, allows one to study whether and when to treat 

■ - - ' this quark as a parton. This reaction therefore sheds light on some of the most 

fundamental aspects of perturbative QCD. We discuss the above issues and review 
the feasibility of their experimental investigation in the light of a large integrated 
luminosity. 
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1 Introduction 



Since the previous HERA workshop in 1991 significant progress has been made on the the- 
oretical side in understanding the production of heavy quarks in electron proton colhsions. 
Improvements in available experimental techniques and particularly the expected increase in 
luminosity amply justify this effort. In general the progress consists of the calculation of all 
0{as) corrections to the processes of interest, thus improving the accuracy of the theoretical 
predictions both in shape and normalization. At the time of the previous workshop the only 
NLO calculations available were for the case of inclusive photoproduction fl]]. In the mean- 
time NLO calculations have also been performed for inclusive electroproduction 0, ^, and 
both have been extended to the fully differential cases |, 0. Therefore, meaningful and 
extensive comparisons between theory and data can now be made. In what follows we review 
how the deeply inelastic electroproduction process allows us to explore, in detail, three areas 
of perturbative QCD in particular. 

We first discuss the inclusive case, via the structure function ^2(0;, Q^, m^). We show that 
this structure function for the case of charm suffers from only very modest theoretical uncer- 
tainty, that its NLO corrections are not too large, and that it is sensitive to the shape of the 
small-x gluon density. Next we treat single particle differential distributions in the charm kine- 
matical variables, and also charm-anticharm correlations. Because many distributions can be 
studied, many QCD tests can be performed. Examples are tests of the production mechanism 
(boson-gluon fusion), studies of gluon radiation patterns, and dependence on scales such as 
deep- inelastic momentum transfer Q, the heavy quark mass m (with enough luminosity one 
can detect a sizable sample of bottom quarks), the transverse momentum of the charm quark, 
etc. Finally, in the last section, we review the theoretical status of the boson-gluon fusion 
description of charm production at small and very large Q. In essence, it involves answering 
the question: when is charm a parton? 



2 Structure Functions and Gluon Density 

This section has some overlap with the more detailed review on heavy flavour structure functions 
in the structure function section. Here we only present the most salient features. The reaction 
under study is 

e-(pe) + Pip) ^ e'ip:) + Qipi)iQipi)) + X , (1) 

where P{p) is a proton with momentum p, Q{pi) {Q{pi)) is a heavy (anti)-quark with momentum 
Pi {Pi = '"^^) s-^d X is any hadronic state allowed. Its cross section may be expressed as 



where 



q = Pe-Pe 



;i + (1 - y)')F2{x, Q\ IT?) - y^FL{x, Q\ ir?) 



X 



y 



p-q 
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2p-q - p-pe 

The inclusive structure functions F2 and F^ were calculated to next-to-leading order (NLO) in 
Ref. 1^. The results can be written as 



Fkix,Q\m') 
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+ E [4 /^^)(cS + cS In 4) + 4, /.(-, /i^)4!] ] ] , (4) 



i=q,q 



where k = 2,L and the upper boundary on the integration is given by z^ax = Q^/{Q^ + 
4m^). The functions fi{x^ii^) , (i = g,q,q) denote the parton densities in the proton and fi 
stands for the mass factorization scale, which has been put equal to the renormalization scale. 
The c2(r7, ^) , c'Hir], > = 5- , 9 , 9 ; 0, 1) and 0> = 5 > ^ ; ^ = 0, 1) are coefficient 

functions and are represented in the MS scheme. They depend on the scaling variables rj and 
^ defined by 

s 

where s is the square of the cm. energy of the virtual photon-parton subprocess which implies 
that in z = Q"^ /{Q"^ + s). In eq. we distinguished between the coefficient functions with 
respect to their origin. The coefficient functions indicated by cl\{V}Oj'^k\{''ljO originate from 
the partonic subprocesses where the virtual photon is coupled to the heavy quark, whereas the 
quantity <ifci(?7,0 comes from the subprocess where the virtual photon interacts with the light 
quark. Hence the former are multiplied by the charge squared of the heavy quark e|^, and 
the latter by the charge squared of the light quark respectively (both in units of e). Terms 
proportional to chGl integrate to zero for the inclusive structure functions. Furthermore we 
have isolated the factorization scale dependent logarithm ln(/x^/m^). A fast program using fits 
to the coefficient functions |^ is available. 

The first thing to note about eq. (Q) is that the lowest order term contains only the gluon 
density. Light quark densities only come in at next order, and this is the reason F2{x,Q'^ ,w?') 
is promising as a gluon probe. To judge its use as such, we must examine some of the charac- 
teristics of this observable. These are: the size of the 0{as) corrections, the scale dependence, 
the mass dependence, its sensitivity to different gluon densities, and the relative size of the 
light quark contribution. These are the issues we investigate in this section. We take the charm 
mass 1.5 GeV, the bottom mass 5 GeV, the factorization scale equal to -^Q'^ + and choose 
at NLO the CTEQ4M [^] set of parton densities, with a two-loop running coupling constant 
for five fiavors and A = 202 MeV, and at LO the corresponding CTEQ4L set, with a one-loop 
running coupling with five fiavors and A = 181 MeV. In Fig. |1| we display F2{x,Q'^ ,m?) vs. x 
for two values of at LO and NLO. The scale dependence is much reduced by including the 
NLO corrections (when varying /i from 2 to 1/2 times the default choice, the structure function 
varies from, at LO, at most 20% and 13% at = 10 and 50 GeV^ respectively, to at most 5% 
and 3% at NLO), but the dominant uncertainty is due to the charm mass and stays roughly 
constant, amounting at NLO maximally to about 16% for = 10 GeV^ and 10% for = 50 
GeV^. The feature that the LO result is mostly larger than the NLO ones is due to the use 
of LO parton densities and one-loop as, and scale choice. Had we used NLO densities and a 
two-loop as, or chosen the scale fi equal to m, the LO result would have been below the NLO 
result. In the first case the size of the corrections is then about 40% at the central values at 

= 10 GeV^, and 25% at = 50 GeV^, and in the second case, at small x, about 20% and 
30% respectively. In the next figure. Fig. 0, we show for the same values of an important 
property, namely the sensitivity of the NLO F2 to different parton density parametrizations. 
In this case we compare the CTEQ2MF set [|T^], whose gluon density stays quite flat when x 
becomes small, and the GRV94 set ||ll[, which has a steeply rising gluon density. One sees 
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Figure 1: F2{x,Q'^ ,m^) vs. x at LO and NLO for two values of Q"^ . The shaded areas indicate 
the uncertainty due to varying the charm mass from 1.3 to 1.7 GeV. 




Figure 2: F2{x,Q^,m'^) vs. x at NLO for two choices of parton densities. The shaded areas 
again indicate the uncertainty due to varying the charm mass from 1.3 to 1.7 GeV. 
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that the difference is visible in the structure function. Finally we remark that the contribution 
of light quarks to the charm structure function is typically less than 5%. The bottom quark 
structure function is suppressed by electric charge and phase space effects and amounts to less 
than 2% (5%) at = 10 (50) GeV^ of the charm structure function. Previous investigations 
of the scale and parton density dependence of F2 using the same NLO computer codes are 



available in [12] and [|T^. 

We conclude that ^2(0;, Q^,m^) for charm production is an excellent probe to infer the 
gluon density in the proton at small x. The NLO theoretical prediction suffers from fairly little 
uncertainty, and the QCD corrections are not too large. See the section on structure functions 
in these proceedings for many more details, where also a comparison with (preliminary) data 
is shown. Therefore in view of a large integrated luminosity, a theoretically well-behaved 
observable, and promising initial experimental studies |T^, |T^] a precise measurement at HERA 
of the gluon density should be possible. 



3 Single Particle Distributions and Heavy Quark Corre- 
lations 

In this section we leave the fully inclusive case and examine in more detail the structure of the 
final state of the reaction 

e-(pe) + P{p) ^ e-(p^) + Q(pi) + g(p2) +X. (6) 

By studying various differential distributions of the heavy quarks we can learn more about the 
dynamics of the production process than from the structure function alone. 

Single particle distributions dF2{x,Q'^,m'^,v)/dv, where v is the transverse momentum pt 
or rapidity y of the charm quark, were presented in NLO in |^ for various choices of x and 
Q^. The LO distributions differed significantly from the NLO ones, so that the effect of 0{as) 
corrections on such distributions cannot be described by a simple K-factor. 

The 0{as) corrections to Fk{x, Q^, m^) in a fully differential form were calculated in Ref. |^ 
using the subtraction method. Recently 0, these fully differential structure functions were 
incorporated in a Monte-Carlo style program resulting in the 0{as) corrections for reaction 
d^). The program for the full cross section, generated according to Eq. (0), allows one to study 
correlations in the lab frame. The phase space integration is done numerically. Therefore, it 
is possible to implement experimental cuts. It furthermore allows the use of a Peterson type 
fragmentation function. For details about the calculational techniques we refer to Ref. |^. 
Here we show mainly results. 

Shown in Fig. ^ are various distributions da/dv for the reaction (^, where the heavy 
(anti)quark has fragmented into a D* meson, with v representing (a) the D* transverse mo- 
mentum p^* (b) its pseudorapidity r]^* (c) the hadronic final state invariant mass W (d) 
for the kinematic range 5 GeV^ < < 100 GeV^, <y < 0.7, 1.3 GeV < p^* < 9GeV and 
\ri^*\ < 1.5. The data are from a recent ZEUS analysis |T5[. The NLO theory curves have 



been produced by using the GRV |11| parton density set, with Peterson fragmentation |16 



The dashed line is for /i = 2m, m = 1.35 GeV and e = 0.035, whereas the solid line is for 



= 2\/Q^ + 4m^, m = 1.65 GeV and e = 0.06. From Fig. H and studies in Ref. [14| it is 
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Figure 3: Differential cross sections and ZEUS data. 



clear that a wide range of studies can be and are being performed already at the single particle 
inclusive level. Preliminary conclusions |1l4| , |15| are that the data follow the shape of the NLO 
predictions quite well, but lie above the theory curves. The HI collaboration |T^ has recently 
shown clearly from the dlna/dxD distribution that the charm production mechanism is indeed 
boson-gluon fusion, (after earlier indications from the EMC collaboration [|T7|) as opposed to 
one where the charm quark is taken from the sea. Here xd = '2\pd*\/W in the 7*P cm. frame 

Next we examine a few charm-anticharm correlations. At the experimental level such cor- 
relations are more difficult to measure since it requires the identification of both heavy quarks 
in the final state. However, with the expected large luminosity that both ZEUS and HI will 
collect, such studies are likely to be done. As an example we show in Fig. |^ the px distribution 
of the pair, p^, and the distribution in their azimuthal angle difference, A(f)'^'^ in the 7*P cm. 
frame for a particular choice of x and Q^. For these figures we used the MRS A' densities |TB . 
Both distributions are a measure of the recoiling gluon jet. 

In summary, differential distributions of deep-inelastic heavy quark production offer a rich 
variety of studies of the QCD production mechanism. Fruitful experimental studies, even with 
low statistics, have been done fl^, |1^ , and with a large integrated luminosity we therefore fully 
expect many more. We finally point out that besides a LO shower Monte Carlo program |jl9|, 
now also a NLO program is available for producing differential distributions. 



4 When is Charm a Parton? 

We return to the inclusive case to ask the fundamental question in the title. The question can 
be more accurately phrased as follows: intuitively one expects that at truly large the charm 
quark should be described as a light quark, i.e. as a constituent parton of the proton, whereas 
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Figure 4: Differential distributions dF2{x , Q"^ , m'^ , p^c) / dp^ and dF2{x , Q"^ , m'^ , Pcc) / dAcj)'^^ at 
X = 0.001 and = 10 GeV^ (solid) and 100 GeV^ (dashed). 

at small (of order m^) the boson-gluon fusion mechanism, in which the charm quark can 
only be excited by a hard scattering, is the correct description. This has been demonstrated 
recently by HI |]14| and ZEUS in [|l^. In this section we examine where the transition between 



the two pictures occurs. 



At LO this issue was investigated in |2y]. A picture that consistently combines both de- 
scriptions, the so-called variable flavor number scheme, is presented and worked out to LO in 
PT| . Here we exhibit where the transition occurs at NLO ||22[| . In other words we will locate the 



onset of the large asymptotic region, where the exact partonic coefficient functions of 0] are 
dominated by large logarithms ln(Q^/m^). These logarithms are controlled by the renormal- 
ization group, and, when resummed, effectively constitute the charm parton density. Here we 
however restrict ourselves to the onset of the asymptotics. Let us be somewhat more precise. 
In m we can rewrite e.g. all terms proportional to as 



X l~ -{n-. f^')H,,{z, % ^) + gC-, f^')H,,{z, % ^)} (7) 



where G{x,fj,'^) is the gluon density and yU^) = '}2i=q^q fi{x, jj?) is the singlet combination 
of quark densities. In the asymptotic regime one may write 

The effort lies in determining the coefficients af'j . Similar expressions hold for the other coef- 
ficients in (1). Taking the limit of the coefficients in is extremely complicated. Rather, a 
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trick |^2[ was used, exploiting the close relationship of the InlQ^/m^) logarithms with coUinear 



(mass) singularities. The ingredients are the massless two-loop coefficient functions of p3[ and 
certain two-loop operator matrix elements. The trick, dubbed "inverse mass factorization". 



X = 10-' 
X = 10--' 
X = 0.01 
X = 0.1 



Figure 5: Ratio of the asymptotic to exact expressions for F2{x, Q^, m?) for the case of charm. 

essentially amounts to reinserting into the IR safe massless coefficient functions the coUinear 
singularities represented by the logarithms ln(Q^/m^). See for details. 



There is another advantage to obtaining the asymptotic expresssions. The terms in eq. (§) 
proportional to e\ have been integrated and full analytical expressions for them exist but 
in the other terms in eq. @) two integrals still need to be done numerically. Therefore in the 
large region the asymptotic formula is able to give the same results much faster, as the 
latter formula needs no numerical integrations. 

In Fig. § we show the ratio of the asymptotic to exact expressions for -F2(x, Q^, m?) for the 
case of charm as a function of for four different x values. Here the GRV |TT| parton density 
set was used, for three light flavors. We see that, surprisingly, already at of order 20-30 GeV^ 
the asymptotic formula is practically identical to the exact result, indicating that at these not 
so large values, and for the inclusive structure function, the charm quark behaves already 
very much like a parton. This is in apparent contradiction with the findings [Q, mentioned 



in the previous section, that the production mechanism is boson-gluon fusion, and illustrates 
that, interestingly, the question in the title can have a different answer for inclusive quantities 
than for differential distributions having multiple scales. 

We finally note that with the results shown in this section also the ffist important step is 
made for extending the variable flavour number scheme to NLO. 
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5 Conclusions 



In the above we have reviewed the many interesting facets of deep-inelastic production of heavy 
quarks. The possibihty of selecting the heavy quarks among the final state particles affords a 
window into the heart of the scattering process, and allows tests and measurements of some 
of the most fundamental aspects of perturbative QCD: the direct determination of the gluon 
density, many and varied studies of the heavy quark production dynamics, and insight into how 
and when a heavy quark becomes a parton. 
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